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Abstract. The high-pressure behaviour of PbS was investigated by angular dispersive X-ray powder diffrac- 
tion up to pressures of 6.8 GPa. Experiments were accompanied by first principles calculations at the density 
functional theory level. By combining both methods reliable data for the elastic properties of rock-salt type 
a- and high-pressure /3-PbS could be obtained. /3-PbS could be determined to crystallise in the CrB-type 
(B33), with space group Cmcm. The reversible ferro-elastic a/ f3 transition is of first order. It is accompa- 
nied by a large volume discontinuity of about 5% and a coexistence region of the two phases. A gliding 
mechanism of {001} bilayers along one of the cubic (110) directions governs the phase transition which can 
be described in terms of group/subgroup relationships via a common subgroup, despite its reconstructive 
character. The quadrupling of the primitive unit cell indicates a wave vector (0,0, ir/a) on the zi-line of 
the Brillouin zone. 

PACS. 91.60.Gf High-pressure behaviour - 61.50.Ks Crystallographic aspects of phase transformations; 
pressure effects - 77.84.Bw Elements, oxides, nitrides, borides, carbides, chalcogenides - 62.20.Dc Elas- 
ticity, elastic constants 



1 Introduction 

Among the IV- VI narrow gap semiconductors lead sul- 
phide, PbS, has been a subject of interest during the last 
three decades, owing to its potential technological rele- 
vance, e.g. in long-wavelength imaging, as diode lasers jl| 
or thermovoltaic energy converters. Many of the experi- 
mental and theoretical investigations have been focused on 
the electronic structure ■ Currently, PbS attracts new 
interest because of its importance for the understanding 
of the complex misfit layer compounds (MX)i +:E (TX2) m 
(M=Sn, Pb, Sb, Bi, rare earth; T=Ti, V, Cr, Nb, Ta; 
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X=S, Se; 0.08<x<0.28; m=l-3) [g, of which PbS can 
be regarded as a component structure. These composite 
structures are single-crystal multi-layer systems, exhibit- 
ing strong anisotropic material properties. 

At ambient conditions a-PbS crystallises in the rock 
salt structure type (Bl) || with the lattice parameter 
a=5. 9240(4) A and space group Fm3m. The elastic con- 
stants of a-PbS were determined by ultrasonic velocity 
measurements ||. Resistivity measurements under high 
pressure on several lead chalcogenides revealed the exis- 
tence of a phase transition Q which occurs in PbS at 2.2 
GPa. First structural studies under high pressure were 
performed using X-ray film methods Later, the com- 
pression behaviour of PbS was studied by means of high- 
pressure energy dispersive X-ray powder diffraction (HP- 
EDX) up to 35 GPa §,|o|Jl|,|lf,|l| . 

The phase transition at 2.2 GPa to a /3-phase having 
orthorhombic symmetry was characterised as a first order 
transition [p| jiqjll| . However, the structure of this phase 
is still a matter of debate. Two different structure types 
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were proposed, viz. the GeS structure type (B16) with 
the space group Pbnm glllUllf and the closely related, 
but higher symmetric CrB structure type (B33) with the 
symmetry Cmcm jlO 11 . A definitive decision between 
these two possibilities has not been possible until now on 
the basis of the experimental data available. 

A second phase transition at about 25 GPa was de- 
tected by HP-EDX and assigned to the transformation 
into cubic 7-PbS with the CsCl (B2) type structure 

Since the equation-of-state parameters for ce-PbS, the 
symmetry of the [3 phase and its compression behaviour 
are hitherto unknown it was found worth reexamining the 
PbS system under pressure. To elucidate the structure of 
the intermediate (3 phase, angle dispersive X-ray powder 
diffraction and density functional theory (DFT) calcula- 
tions were performed. Here the results of the study of the 
crystal structure of PbS under high pressure are presented. 
The mechanism of the phase transition is discussed using 
the symmetry relations between the space groups of the 
two phases. 



2 Experimental methods 

Crystals were grown by chemical vapour transport (CVT) 
in evacuated quartz tubes in a gradient furnace (1053 to 
923 K) with stoichiometrically mixed elements as start- 
ing material and iodine as transport agent. In the low- 
temperature zone cube shaped crystals were obtained af- 
ter several days with edge lengths up to 3 mm. The chem- 
ical composition of the crystals was checked by electron 
micro-probe analysis and the phase purity was studied by 
conventional X-ray powder diffraction. 

The high-pressure experiments were performed us- 
ing diamond anvil cells Ej]. NaCl was used as internal 
standard for the pressure determination by the Decker 
equation-of-state Cq]. Several runs were performed us- 
ing a 4:1 methanolTethanol mixture for maintaining hy- 
drostatic conditions. The highest pressure obtained was 
6.8 GPa. Angular dispersive diffraction data were collected 
with Mo- if a radiation on a Mar2000 image plate diffrac- 
tomcter. The geometry correction for the radial integra- 
tion of the two-dimensional data and the transformation 
into standard one-dimensional powder patterns were per- 
formed using Fit2d [lql Whole powder pattern fitting 
employing Fullprof |fL7| was performed for the precise 
determination of the unit cell parameters. 



3 Computational aspects 

Geometry optimisations were performed for a-PbS and 
for both proposed structur e ty pes of /?-PbS. The ab initio 
plane- wave CASTEP code pj,hj was used which is based 
on the density functional theory (DFT). The exchange 
correlation interaction was accounted for by using the 
Perdew, Burke and Ernzerhof version of the generalised 
gradient approximation (PBE-GGA) of the exchange- 
correlation functional J2^,||l[. Ultrasoft pseudopotentials 




Fig. 1. Pressure dependence of diffraction patterns of PbS 
from room pressure to 6.87 GPa. The patterns are shifted with 
increasing pressure. Miller indices are given for Q-PbS (lower 
row), NaCl (middle) and, /3-PbS (top row). 



[|22j|23|] with a maximum cutoff energy of the plane waves of 
290 eV were used to describe the electron-ion interaction. 
The integration over the Brillouin zone was performed 
employing the Monkhorst-Pack |24| sampling scheme of 
reciprocal space, with a distance of 0.04 A -1 between 
the sampling points. All structural parameters not con- 
strained by the space group symmetry were relaxed. 

The calculation of the elastic constants is based on 
applying small strains to the already relaxed structure, 
followed by relaxation of the atomic positions. In order to 
avoid anharmonic effects the maximum applied distortion 
was about 0.2 %. Several amplitudes for each strain com- 
ponent were used. Elastic constants were determined from 
a linear fit of stress to the applied strain. 



4 Results and discussion 

Some integrated powder patterns are presented in figure 
|l]. The onset of a phase transition was detected at 2.5 GPa 
by the appearance of additional reflections. The new phase 
could be indexed with an orthorhombic unit cell and cell 
parameters a=3.98(2) A, 6=11.11(7) A and c=4.16(4) A. 
Coexistence of the two phases was observed up to a pres- 
sure of 6 GPa. 

The powder patterns obtained above 2.5 GPa could be 
indexed with the same unit cell for both, B16 and B33, 
structure types. A structure refinement of the /3-phase in 
either case turned out to be impossible because of the 
strong peak overlap with reflections stemming from the 
gasket material (Inconel steel). Thus, a discrimination be- 
tween both possible structure types was not possible by 
the diffraction method. 



4.1 Lattice parameters and elastic properties of a- 
and /3-PbS 

The cell parameters of the orthorhombic /3-phase are re- 
lated to those of the cubic a-phase by a = \/2a c /2, b = 
2a c , c — y/2a c /2, with subscripts o and c referring to 
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Table 1. Parameters of Birch Murnaghan equations-of-state 
derived from experimental and calculated data and values 
of the linear compressibility along the directions of the lattice 
vectors for ct- and /3-PbS. The linear compressibilities from the 
DFT calculations were derived from the elastic compliances as 
k a = sn + S12 for the cubic system and k a = su + s\2 + S13, 
kb = 322 + S12 + S23, k c — S33 + S13 + S23 for the orthorhombic 
system ^q| . 





exp. 


calc. 


b [GPa] 


51.0(1.2) 


52.43(9) 


a-PbS b' 


4.3(9) 


4.69(2) 


k a [GPa" 1 ] 


0.0063(8) 


0.00628 


b [GPa] 


30.9(4) 


25.6(6) 


b' 


4(fixed) 


6.9(2) 


/3-PbS k a [GPa" 1 ] 


0.0049(1) 


0.0069 


k b [GPa" 1 ] 


0.0032(4) 


0.0095 


k c [GPa- 1 ] 


0.0055(2) 


0.0024 



the orthorhombic and the conventional cubic lattices, re- 
spectively. Figure ^a shows the dependence of (pseudo-) 
cubic lattice parameters on pressure for both phases. 

The normalised volume changes are given in figure |^b. 
Because of the narrow pressure range observed in the ex- 
periment the volume at zero pressure Vo of /3-PbS could 
not be obtained unambiguously from the equation-of-state 
(eos) fit. Instead, the volume at zero pressure from the ab 
initio calculations was used to normalise the experimental 
data of the /3-phase. 

The non-zero components of the calculated elastic stiff- 
ness tensor (cy-) (in Voigt's notation, units GPa) are 
en = 141(3), C12 = 34(2) and c 44 = 19.97(5) for a-PbS at 
ambient pressure. They agree well with the available ex- 
perimental data (en = 149, c\i — 35 and C44 = 29 GPa) 
determined at 5 K Q. For /3-PbS no experimental data 
are available, the computed values at 4 GPa are 

/ 105(5) 11(6) 67(5) \ 
89(8) 30(3) 
103(7) 

21(1) 

71.02(6) 

V 12(2)7 

Linear bulk moduli determined from elastic constants are 
70(2) GPa for a-PbS and 53(3) GPa for /3-PbS. The value 
of the linear volume bulk modulus of the cv-phase is higher 
than that of the high-pressure /3-phase. This reflects the 
diffe ren t compression behaviour in the two phases (sec- 
tion (L3) . The bulk moduli bo and their pressure dependen- 
cies V were determined for a- and /3-PbS by fitting Birch 
Murnaghan equations-of-state |25| to the normalised vol- 
ume data. The fit results are given in table [l] together with 
the linear compressibility data k. 

The bulk modulus bo, as well as the pressure depen- 
dence of the bulk modulus b', for a-PbS show excellent 
agreement between experimental and calculated values. 
Also, the linear compressibilities derived from a linear 
fit to experimental data and from calculated compliances 
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Fig. 2. (a) Pressure dependence of the experimentally deter- 
mined cell parameters for a- (□ : a c ) and /3-PbS ( A : V2c , 
o : \/2a a , o : b /2 ). (b) Pressure dependence of the normalised 
unit cell volume of a-(°) and /3-PbS (<>). The lines in (a) cor- 
respond to linear fits to the experimental data points and in 
(b) to fits of a 3 rd -order Birch-Murnaghan equation-of-state to 
the ab initio data. 

agree excellently. The limited pressure range of the exper- 
imental data prevented us from determining b' for /3-PbS. 
Hence, the bulk modulus was determined with a second 
order Birch Murnaghan equation-of-state (b 1 :— 4). It is 
well known that bo and b' are highly correlated variables 
in eos-fitting Mj. Consequently, bo is higher in the ex- 
periment, since b' is smaller comparing the numbers be- 
tween experiment and theory. The resulting bulk modulus 
is smaller than data reported recently [ |ll| . However, di- 
rect comparison is difficult since the authors of |ll[] used 
a different equation of state. 



4.2 Stability of the (3 phase 

As stated earlier, the quality of our diffraction data did 
not allow the unambiguous assignment of the symmetry 
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Fig. 4. (a) Variation of the Pb-S bond lengths in /3-PbS from the ab initio calculations. In (b) and (c) coordination polyhedra 
and the crystal structure of /3-PbS are drawn. Applying the labels given for the S atoms of the single-capped trigonal prism, 
shown in (c) the following bonds can be distinguished in graph (a): squares represent the Pb-S(l) bond pointing at the apex 
of the pyramid located within the bilayer in (b), circles correspond to Pb-S(2-5) bonds which constitute the base plane of the 
pyramid and one face of the prism, and triangles denote the Pb-S(6,7) bonds between adjacent bilayers and, being located 
inside the prism. 



and, thus, discrimination between both possible structure 
types B16 and B33 for /3-PbS. By contrast, the quantum 
mechanical calculations made the distinction possible, be- 
cause only for the CrB-type structure (B33) convergence 
was achieved. The enthalpy for this phase with respect 
to the a-phase is shown in figure |as a function of pres- 
sure. A transition pressure of 9.6(1) GPa for the a/ (3- 
transformation follows from the crossover of the fitted 
curves. This transition pressure is higher than the experi- 
mentally determined values (onset of the transition at 2.6 
GPa, this study; 2.2 GPa §). However, DFT-GGA cal- 
culations such as those presented here suffer from system- 
atic 'underbinding' and the neglect of temperature effects. 
Furthermore, they refer to perfect ideal crystals. Hence, a 
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Fig. 3. Pressure dependence of the molar enthalpy difference 
between a- and (3- PbS. 



discrepancy of 7 GPa between observed and calculated 
transition pressures can be tolerated. Also, the main aim 
of such calculations is not to reproduce numbers, but in- 
stead provide insight into the transition mechanism, and 
for this purpose such calculations are well suited. 



4.3 Compression mechanism 

In ev-PbS the lead atoms are sixfold coordinated by sul- 
phur atoms in a regular octahedral arrangement. The only 
structural parameter in the Bl structure type not con- 
strained by the space group symmetry is the magnitude 
of the cell parameter. The Pb-S bond length is a c /2 and, 
of course, the compression mechanism results in shorten- 
ing of the bond length. 

No major change in the linear compressibility be- 
tween cubic and orthorhombic PbS follows from the ex- 
perimentally determined cell parameters (Fig. ^d). How- 
ever, Pb-S bond lengths derived from the ab initio cal- 
culations (Fig. |^a) clearly show that the compression 
behaviour is anisotropic. Layered chalcogenides show a 
strongly anisotropic compression with the strongest com- 
pression perpendicular to the layers (e.g. TiS2 ]2q] , SnS2 
[ p9[ , NbS2 ppf ). The main compression mechanism in 
these compounds is narrowing of the van-der-Waals gaps, 
whereas the layers themselves are bonded predominately 
ionic. The structure of /3-PbS can be thought of as a 
layer structure as well, consisting of rock salt like bilay- 
ers, stacked along the 6-axis of the orthorhombic unit cell. 
This is shown in figure ^b where two simple coordination 
polyhedra can be distinguished. Within the bilayers Pb is 
found in a fivefold coordination for which a square pyra- 
mid is an idealised geometry. The basal pinacoid face of 
the square pyramid forms together with two sulphur atoms 



K. Knorr et al.: The high-pressure a/f3 phase transition in lead sulphide (PbS) 



5 



Table 2. Geometric conditions for an ideal single-capped trig- 
onal prism Five different edges with index i have to be 
distinguished occurring with frequency m and different edge 
lengths di. The radius of the circumsphere is r and labels are 
given according to Figure fy. 



i 


di 


" i 


label 


1 


1.277r 


4 


Sl-S2,3,4,5 


2 


1.297r 


2 


S2-S3, S4-S5 


3 


1.233r 


4 


S2-S6, S5-S6, S3-S7, S4-S7 


4 


1.195r 


1 


S6-S7 


5 


1.478r 


2 


S2-S5, S3-S4 



Table 3. Distortion measures of the single-capped trigonal 
prism in /3-PbS at selected pressures, obtained from DFT cal- 
culations. 



p [GPa] 


Sy [A] 


f w [A] 


ADM 


A d [A] 





.3594 


3.201 


.0289 


.0702 


2 


.2149 


3.117 


.016 


.0512 


7 


.1005 


3.019 


.0114 


.0416 


10 


.0692 


2.979 


.0111 


.0405 


15 


.0355 


2.928 


.0102 


.0382 


21 


.0095 


2.879 


.0105 


.0381 


25 


.0099 


2.852 


.0099 


.0367 



from the next bilayer an almost trigonal prism. Figure |]a 
shows the strongest compression for this particular co- 
ordination polyhedron. The Pb-S(6,7) bonds connect the 
layers. The two different bonds within the bilayers are less 
compressible. While the pyramids deform almost isotrop- 
ically, the main compression mechanism is a distortion of 
the trigonal prism, due to a change of the Pb-S(6,7) bond 
lengths. This mechanism effectively influences the length 
of the &-axis. Changes in the a and c axes are related to 
lateral distortions within the bilayers, i.e. changes of the 
Pb-S(2-5) bond lengths, and of the angles between them. 

Admittedly, a description of the coordination in /3-PbS 
by two different polyhedra is not fully appropriate since 
this approach doesn't account for the increase of the co- 
ordination number with pressure. Alternatively, the coor- 
dination of the lead atoms can be described as sevenfold 
with the geometry of a single-capped trigonal prism fl3j| . 
The capped trigonal prism in /3-PbS is shown in Figure^jc. 
The ideal geometry of a single-capped trigonal prism re- 
quires that all apices are located on the surface of a sphere 
and the edges di fulfil the conditions given in table |^ J3lj . 
The deviation of this coordination polyhedron from ideal 
can be quantified by applying e.g. the Pinsky and Avnir 
continuous-symmetry-measure (CSM) |3^| . Here a simpli- 
fied CSM formulation was used. 

From the five different edges di radii r, were deter- 
mined according to table || and from the weighted average 
radius 

5 5 

1=1 2 = 1 

the edge lengths of the corresponding ideal polyhedron d* 
could be calculated. The average distortion measure ADM 
is then given by 



ADM 



Eli Mf 



100 



(2) 



and the mean edge length distortion is 



A d = 



\ 



j> (di-d*f/Yl 



m. 



(3) 



i=l 



The ADM is normalised with respect to the edge lengths of 
the ideal polyhedron in order to compensate for the effect 



of the compression. Hence, it is a dimensionless number. 
The compression of the coordination polyhedron is repre- 
sented by the change of of the circumsphere's radius f w 
with pressure. The mean edge length distortion Ad is a 
measure of the variance of the edges. The resulting values 
for selected pressures are given in Table ^. At low pressures 
the conformation of the coordination polyhedron deviates 
most strongly from that of an ideal single-capped trigo- 
nal prism |31j (ADM=0). With increasing pressure ADM 
and Ad indicate increasing regularity of the coordination 
polyhedron. Above approximately 10 GPa some residual 
deformation remains virtually constant. The position of 
the central lead atom does not coincide with the conjoint 
centroid of the seven sulphur ligands (the difference 8 y 
is along the orthorhombic b direction and given in table 
||) . With increasing pressure the position of Pb tends to- 
wards that of the sulphur's centroid. A similar behaviour 
has recently been reported for a-PbF2 |||, where Pb is 
elevenfold coordinated by F in a tri-capped trigonal prism. 



4.4 Characterisation of the a/f3 phase transition 

Coexistence of the two phases between 2.5 and 6 GPa was 
observed. Since the peak profile remained sharp through- 
out the whole pressure range covered by the experiment, 
pressure inhomogeneities can be excluded as the origin of 
this phase coexistence. Instead, it reflects the first order 
character of the transition (e.g. |]l0|Jl if] ) . This is further 
supported by the observed discontinuity in the molar vol- 
ume of 4.7 % at the phase transition. 

The transition is reconstructive. This is justified by i) 
the change of the coordination number from six to seven, 
and ii.) because there is no direct group/subgroup rela- 
tionship between the two space groups Fm3m and Cmcm. 
Figure p3 shows the symmetry relations in terms of klassen- 
gleiche (index k) and translationengleiche (index t) sub- 
groups together with the changes of the unit cell dimen- 
sions. The crucial step here is the transition to space group 
Pbnm realised by a displacement < Aa < 1/2 of every 
second (OOl)-bilayer along the orthorhombic a-, i.e. cubic 
(110), direction. The space group Cmcm is obtained as a 
limit state, when Aa takes the value of 1/2. It is interest- 
ing to note that the displacement is the same as for the 
preferential glide system (110){001} of NaCl. In the (llO) 
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FmSm 



a-PbS 



t.l 



(a-b)/2, (a+b)/2, c 
14/m.nim 

tg 
Immm 



the MX subsystem in misfit layer compounds. Of course, 
for a better understanding of the elastic properties of the 
MX-TX2 multi-layer systems further investigations are re- 
quired. 
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Fig. 5. Group/subgroup diagram for the a//3 phase transi- 
tion in PbS, showing a possible transition pathway from space 
group Fm3m to Cmcm via the common subgroup Pbnm. 



directions easy displacements of {001} layers are possi- 
ble because movement of the ions occurs without having 
contact of equally charged ions. 

The a/P transition in PbS is accompanied by a qua- 
drupling of the unit cell. The wave vector of the transition 
is (0, 0, it /a) in the midpoint of the Zi-line connecting the 
-T-point at the zone centre with the X-point on the bound- 
ary of the Brillouin zone of the cubic face-centred lattice. 



5 Summary and conclusion 

The elastic behaviour of a- and /3-PbS was measured up to 
6.8 GPa. Bulk moduli, elastic constants and linear com- 
pressibilities were determined by DFT calculations and 
compared to experiment. The symmetry of /3-PbS was 
found to be Cmcm. The differences in the compression 
behaviour between a and /3-PbS can be understood on 
the basis of the different crystal structures and conse- 
quently, different compression mechanisms. PbS bilayers 
as found in the /3-phase occur as MX subsystem of the 
chalcogenide misfit layer compounds [ Ij . Structural inves- 
tigations under high pressure are not known for the misfit 
layer compounds. Therefore, our findings for the compres- 
sion mechanism in /3-PbS may serve at least as model for 
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